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Abstract. The paper presents a study of two full-core, edge dislocations of opposite 
Burgers vectors in 4H-SiC, conducted using the first-principles density functional 
theory methods. We have determined the creation energy of the dislocations as a 
function of distance between their cores. The radial distribution function has been 
applied to examine strong impact of the dislocations on the local crystal structure. 
The analysis of the electronic structure reveals mid-gap levels induced by broken 
atomic bonds in the dislocation core. The maps of charge distribution and electrostatic 
potential have been calculated and the significant decrease of the electrostatic barriers 
in the vicinity of the dislocation cores has been quantified. The obtained results have 
been discussed in the light of available experimental data. 


PACS numbers: 71.30.-fh, 71.38.-k, 64.70.Kb, 75.50.Gg 
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Introduction 

The Silicon Carbide (SiC), a wide band gap semiconductor, is a promising material for 
high-voltage and high-frequency nanoelectronic devices IB El- Very good operational 
quality of SiC results from high values of breakdown voltage (~ 10® V/cm), high 
charge carrier mobility, high temperature stability and high thermal conductivity j3]. 
Additionally, this material has very good mechanical properties and resistance to 
radiation damage. Unfortunately, the electronic properties of epitaxial layers strongly 
depend on the material quality. The presence of intrinsic defects and impurities which 
arise during crystal growth process substantially limit applications of SiC. 

Dislocations are the main crystal defects in SiC. They deteriorate the performance of 
high electric field devices such as Shottky and p-n diodes HE]. Two typical components 
of dislocations with the direction along [0 0 01] are the screw and edge dislocations with 
Burgers vectors [0 0 01] and | [2 11 0], respectively. The standard density of dislocations 
observed in the good quality 4H-SiC epilayers reaches 10^ cm“^ jH]. The electrical 
characteristics of 4H-SiC photodiodes reveal that the screw and edge dislocations reduce 
the breakdown voltage by 3.5% and 2%, respectively, and increase the leakage currents 
compared to systems without crystal defects [7]. Dislocations also influence transport 
properties of 4H-SiC by increasing recombination activity [8] and reducing the diffusion 
length in the material [9|. 

Various mechanisms are responsible for modihcations of electronic structure in the 
presence of dislocations. In the dislocation core, the broken bonds give rise to acceptor 
levels within the band gap and make this region electrically active (TOlllI]. Dislocations 
also act as trapping centres and sources for point defects [12]. In the vicinity of the 
core, the strain field may induce additional states located near the conduction band 
edge [13]. In epitaxial SiC, several gap levels have been detected using the deep-level 
transient spectroscopy (DLTS) and their connection with the intrinsic defects have been 
analysed [HUS]. 

The structural properties of dislocations have been previously studied theoretically 
in very few materials using both empirical potentials [16] and ab initio methods based 
on the density functional theory (DFT) [H EHl IH 120] • The DFT provides tools that 
allow one to study changes in electronic structure induced by dislocations [2T]. Such 
computations revealed that the dislocation core in GaN induces deep-gap states [22l 123] , 
which strongly influence the electronic and optical properties of this semiconductor [21] . 
The point defects such as vacancies and oxygen dopants can be easily trapped at the 
core of the edge dislocations providing additional mid-gap states [221 [21] • The influence 
of partial dislocations on the structural and electronic properties have been investigated 
in the pure SiC [21] and crystals doped with impurity atoms [27]. However, to our best 
knowledge, electronic properties of the full-edge dislocation has not been studied with 
ab initio methods yet. 

In this work we study the changes in the crystal structure and electronic properties 
of 4H-SiC induced by the edge dislocation. We analyse in detail the local lattice 
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distortion by means of the radial distribution function (RDF) - the density of atoms in 
the spherical shell around given atom averaged over whole structure [28l |29l |30] . The 
calculations disclose signihcant redistribution of charges and electrostatic potential in 
the region of the dislocation core. The electronic states located in the insulating gap, 
arising from the dislocation core atoms, show a very weak dispersion in the perfect 
correlation with the distribution of charges in the distorted region. Estimated minimal 
electrostatic barriers for ideal and distorted systems throw a new light on the decrease 
of the breakdown voltage in the defected 4H-SiC crystal. 

Methods 

All calculations have been executed with the Vienna Ab-Initio Simulation Package 
(VASP) [311 132] using the full-potential projector-augmented wave method [331 El] 
within the generalised-gradient approximation (GGA) [331136] EZ]- The following valence 
base conhgurations have been included: Si and G 2s‘^2p^. The integration over 

the k-point space has been performed over the 2x2x2 Monkhorst-Pack mesh [3E] and 
the energy cut-off for the plane waves expansion was equal to 500 eV. The structures of 
ideal and defected crystals have been fully optimised (with respect to lattice parameters, 
stresses, and internal degrees of freedom) in the supercells consisting of 8x6x 1 primitive 
unit cells, containing 384 and 346 atoms, which correspond to 25Axl9AxlOA and 
22Ax19Ax10A crystal volumes, respectively. We have used the conjugate gradient 
technique with the energy convergence criteria set at 10“^ eV and 10“^ eV for the 
electronic and ionic iterations, respectively. 

Results 

Crystal structure 

The pure 4H-SiG crystal structure has been reproduced correctly with the lattice 
constants a = 3.095 A and c = 10.132 A, which are very close to experimental values 
that amount to 3.073 A and 10.053 A, respectively. To fulhl the periodic boundary 
conditions we consider two full-core edge dislocations with opposite Burgers vectors 
I [2 11 0] and |[2 110], inserted into the system using Visual Editor of Grystal Defects 
(VEGDs) [33] • It involves removal of 38 selected atoms (19 SiG molecules) changing 
the total number of them, which has been taken into account in all comparisons to the 
ideal structure. Nevertheless, the relaxed crystal volume noticeably increases (by about 
1%) due to the local destruction of hexagonal close-packed structure. However, since 
the supercell is big enough, the elongation of the average near neighbour bond is below 
0.5% in comparison to its length in the ideal system. 

The creation energy per unit length of a single edge dislocation has been calculated 
as a function of distance between their cores using the following formula [23] : 

e = {E-Eq)/2c, 
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Figure 1. OfF-bottom (upper panel) and top (lower panel) view of the optimised 
crystal structure of 4H-SiC with two full-core edge dislocations. Both dislocation 
edges (I and II) have been distinguished by colour. 


where E is the energy of the system with a pair of dislocations (346 atoms), Eq is 
the gronnd state energy of the undisturbed system normalised to the same number of 
formula units. For relative distance between cores equal to 6.73 A we have obtained 
e = 2.17 eV/A per supercell. Additionally, it has been found that e increases for larger 
distances due to attractive force between two dislocations of opposite sign [40]. The 
Peierls energy barrier between two states of dislocation core at adjacent lattice nodes 
along the common glide plane is very low. Therefore, to prevent the annihilation of 
defects during optimisation process we have trapped both dislocations by pinning their 
cores at different glide planes. 

The optimised 4H-SiC structure with two edge dislocations is presented in Fig. [Tj 
The positions of dislocation lines have been marked by atoms distinguished by different 
colour. Carrying out crystal structure relaxation, without any symmetry constrains, we 
have found that it tends to self-recover - reconstructing as much as possible of pure 
4H-SiC coordinations and so gathering most of the disorder to one of the dislocation 
cores (I in Fig. [^. 

To analyse the influence of the dislocation on its nearest neighbourhood, the 
RDF for the dislocation dipole has been calculated and analysed. In Fig. one can 
immediately notice that the introduction of the edge dislocation disturbs hexagonal 
close-packed structure which results in slight increase of the near neighbour and next 
near neighbour distances visible as shift from pure 4H-SiC positions (dashed lines). We 
discovered that each coordination shell can be fitted with a two-component function 
(two Gaussian lines with parameters collected in Table [^: the central part with small 
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Figure 2. RDF of the first-(left panel) and the second-coordination shell (right panel) 
for the edge dislocation dipole in 4H-SiC (points), fitted with the two-component 
Gaussian (solid lines). The dashed lines denote positions of the coordination shells for 
the ideal 4H-SiC system. 


Table 1. The parameters of Gaussian functions (positions and FWHM) and calculated 
peak areas of two-component fit of RDF presented in Fig.[^ u and d parts corresponds 
to the undistorted and distorted regions, respectively. 



near neighbour 

next near neighbour 

u 

d 

u 

d 

position 

1.911 

1.896 

3.115 

3.121 

area 

0.321 

0.387 

0.154 

0.560 

FWHM 

0.034 

0.101 

0.024 

0.185 


width and the rest reproducing broad tails behaviour. In particular this is clearly visible 
for the second coordination shell (right panel in Fig. |^. 

We have found that the components with a narrow and broad full width at half 
maximum (FWHM) correspond to the undistorted part of the crystal and disordered 
dislocation core, respectively. For bigger supercells with smaller densities of dislocations 
the latter component will remain mostly unchanged while the former one will increase 
and converge towards a delta function. 

Band structure 

Next, we analyse the influence of edge dislocations on the band structure of the 4H-SiC. 
The general observation is that the positions of the main bands in the electron density 
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Figure 3. The electronic density of states (eDOS) calculated for pure 4H-SiC (filled 
contours) compared with results obtained for the distorted system (thin colour lines). 
Fermi energy level has been marked with dashed line. 

of states (eDOS) remain mostly unchanged (see Fig. [^. The characteristic peaks in the 
valence and conduction bands of the ideal structure become broadened and smeared out 
in the system with dislocations. The strongest changes are found inside the energy gap 
of the pure 4H-SiC crystajj} where additional localised states appear. One can hnd sharp 
deep states delivered by atoms with broken nearest neighbours bonds and broadened 
shallow bands close to the top of the valence band and the bottom of the conduction 
band. The presence of shallow bands effectively reduces the magnitude of the insulating 
gap. The states located in the range ~ 0.2 — 0.6 eV below the conduction band have 
also been found in the SiC systems with the stacking faults imsg and Shockley partial 
dislocations |2S]. The latter induce also the occupied states about 0.4 eV above the 
valence band. All deep states are very narrow which indicates their weak dispersion 
in the Brillouin zone and localised character. It is possible to trace the distribution of 
individual states in the whole spectrum, as well as to identify atoms which give their 
contribution to the eDOS in a particular energy range. We have found that the deep 
states originate from the atoms in the distorted region of the dislocation core. The 
detailed analysis has led to the conclusion that also partial eDOS projected to the spd 
states has been modihed. For the conduction band, where in the perfect 4H-SiC d states 
dominate, as well as for the localised states inside the gap, the main contribution for 
distorted system stems from p states. 

f The theoretical gap Eg = 2.3 eV is smaller than the measured value (3.2 eV) due to well-known gap 
underestimation problem in the LDA/GGA approaches. 
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Figure 4. (a) Off-bottom view of charge distribution in distorted system. The slice 
plane crosses through both edge dislocations cores, (b) Top view of electrostatic 
potential derived from charge distribution averaged along c direction. (c) The 
electrostatic potential along c-axis direction in PI (dislocation core), P2 (between 
atoms in the ordered region), and P3 (atomic chain) places marked in panel (b). 
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Electric charge distribution 

The charge density is a well-defined, primary qnantity in the DFT. Fig. presents the 
charge redistribntion aronnd edge dislocations. In particnlar, one can see the elongated 
strnctnres with a very low charge density (ontlined with dotted line in Fig. |^) in the 
slice plane crossing both dislocation cores. On the other hand, the neighbonring atoms 
visible in Fig. have increased density of carriers connected with broken interatomic 
bonds. The effect is more prononnced for the dislocation I, marked with the thick 
solid line in the fignre. The electrostatic potential derived from the charge distribntion 
and averaged over c direction is presented in the top view projection (Fig. |^). The 
decreased charge density regions exhibit lowered energy barriers which allow for easier 
flow of the carriers. Again, the effect is mnch stronger for the dislocation I with the 
bigger disorder aronnd its edge. Additionally, in panel (c) the electrostatic potential 
along the c direction in three different locations of distorted system has been shown. 
Its periodicity reflects ABCB stacking in the 4H-SiC strnctnre. In accordance with the 
intnition, electric charges enconnter the highest barriers along the atomic chain (cnrve 
P3 in Fig. |^) and the lowest throngh the dislocation core (cnrve PI). 

To qnantify the infinence of dislocations on the energy barriers for the electron 
transport along c-axis, we have analysed the local electrostatic potential in the crystal. 
We have taken into acconnt that the path over minimal barriers can wander not 
necessarily along a straight line parallel to the c direction. Therefore, for each pair 
of positions: tq = {x^y^Zo) and ri = {x^y^Zi) - where Zq is at the bottom of the 
periodic nnit while zi is at its top (see Fig. |^) - we have derived a minimal energy path 
and the barrier height for this path from rg to ri. The calcnlation nsed the basin-filling 
segmentation algorithm uni with bisection search for minimal limiting energy level. 
The resnlts indicate snbstantial lowering of the barrier in the vicinity of the dislocation 
cores denoted by PI in Fig. |^. The calcnlated difference between the minimal barrier 
height in the non-distorted area and in the vicinity of the dislocation core is 0.8 V. This 
valne is mnch lower than the differences in barrier heights presented in Fig. dne to 
the nonlinear character of the minimal barrier path. The calcnlated energy barriers are 
obtained with zero voltage bias and as snch are not snfficient to determine qnantitatively 
the charge transport properties of the material bnt are an important step towards fnll 
nnderstanding of charge carriers’ behavionr. 

On the basis of the presented resnlts, the following mechanism of the insnlating 
properties weakening and the breakdown voltage decreasing in the 4II-SiC monocrystal 
may be proposed. Very low energy barriers along glide planes cannot oppose the 
formation of dislocations aronnd point defects and local lattice stresses as well as their 
migration throngh the crystal. In a dislocation core, part of interatomic bonds are 
broken cansing the creation of deep states inside the semicondnctor gap. Additionally, 
crystal distortions modify the atomic potential enforcing the shallow states formation 
and conseqnently narrowing the forbidden gap. Finally, the elongated regions of the 
reduced charge density with fiat electrostatic potential are formed along a dislocation 
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core. The lack of barriers inside so created tunnels, in connection with additional 
states in the energy gap, may enhance a carrier flow through the distorted areas, 
significantly influencing decrease in the breakdown voltage. This scenario provides 
the plausible explanation of the current-voltage characteristics of the 4H-SiC avalanche 
photodiodes, which show a pronounced decrease of the breakdown voltage and increase 
of leakage currents due to a single edge dislocation [7j. Our analysis of the mid-gap 
levels also enables a better understanding of a pronounced impact of dislocations on the 
carrier lifetime [8] and diffusion length [9] , and can be helpful in interpreting the DLTS 
measurements [I1[I5]. Moreover, similarly to other strongly defected systems [H], the 
changes in the electronic structure may be related to observable effects in the lattice 
dynamical and optical properties of SiC [45] . 

Conclusions 

Summarising, in this work we have successfully modelled a pair of edge dislocations using 
ab initio methods. We have confirmed experimental findings concerning small but well- 
visible lattice constants elongation, energy gap narrowing, and electrostatic barriers 
reduction. Furthermore, with atomic-scale resolution, we have explained foregoing 
processes. We have shown that (i) the crystal structure is strongly disturbed in the small 
vicinity of the dislocation core, (ii) additional energy levels occurring in the energy gap 
belong to the atoms with broken bonds occupying the core neighbourhood, (iii) existence 
of spatial tunnels, with atoms delivering localised states to the band structure on its 
sides, signihcantly decreases electrostatic barriers and should be considered as one of the 
primary factors responsible for experimentally observed reduction of breakdown voltage. 
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